Spinon-Holon Attraction in the Supersymmetric t — J Model with l/r 2 -Interaction 
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We derive the coordinate representation of the one-spinon one-holon wavefunction for the super- 
symmetric t — J model with l/r 2 -interaction. This result allows us to show that spinon and holon 
attract each other at short distance. The attraction gets stronger as the size of the system is in- 
creased and, in the thermodynamic limit, it is responsible for the square root singularity in the hole 
spectral function jjj]. 
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The low- lying excitations of one-dimensional (1-d) 
strongly-correlated electron systems are not Landau's 
quasiparticles or quasiholes carrying both charge and 
spin Q, but rather collective modes carrying spin-1/2, 
but no charge (spinons), or charge 1, but no spin (holons) 
Spinons and holons are semions, i.e. particles with 
statistics half that of regular fermions ||,|5| . In the ther- 
modynamic limit, the physics of 1-d correlated electron 
systems is described by a gas of spinons and holons 
The corresponding energy is additive, being the sum of 
the energies of each isolated particle. Although this im- 
plies that the spinon-holon interaction energy is irrele- 
vant in the thermodynamic limit, it gives no information 
about short-distance dynamics. 

In this letter we carefully analyze the short-distance 
interaction between a spinon and a holon in an exact so- 
lution of the supersymmetric t — J model with 1 /r 2 inter- 
action ("Kuramoto- Yokohama" (KY)-modeF) 0. The 
KY-model is a system of electrons on a lattice with pe- 
riodic boundary conditions. Double occupancy of a site 
is forbidden by strong Coulomb repulsion. Unoccupied 
sites are allowed and therefore holes can live on the lat- 
tice. Charge hopping, Coulomb, and spin-spin antiferro- 
magnetic terms are inversely proportional to the square 
of the chord between the corresponding sites. The Hamil- 
tonian is given by 
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where z a = exp(2Tria/N) (a is a lattice site, N is the 
number of sites), S a is the spin operator at site 
the electron operator at site a, n a<J = c^Cac, and P is 
the Gutzwiller projector that annihilates configurations 
with doubly-occupied sites: P — I1q(1 — n a^n a i). At 



filling- 1/2, the KY model coincides with the Haldanc- 
Shastry (HS) model of 1-d antiferromagnet. 
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FIG. 1. Square of the spinon-holon wavefunction |pmn(z)| 
defined by Eq.Q for m = N/2 - 1 and N = 600. The 
probability peaks up at short separation between spinon and 
holon, while it does not depend on the distance at large sep- 
arations. The inset shows the function around the origin 
for N — 200, 400, 600. The value at the origin diverges as 
N — > oo. 

In recent papers ||||, we have worked out the real- 
space coordinate representation for two-spinon wavefunc- 
tions of the HS model. In this letter, we extend our 
technique to define and work out the coordinate rep- 
resentation for one-spinon one-holon eigenfunctions and 
the corresponding Schrodinger equation for i?KY close to 
filling- 1/2, where antiholon excitations, made out of a 
hole created in a holon sea, are ruled out ||,[l(J . Spinon- 
holon interaction and its nature follow from the behav- 
ior of the exact solution of this equation. In Fig|l] we 
plot the result. The probability does not depend on the 
spinon-holon separation when they are far apart. This 
shows that the two particles are noninteracting at large 
distances. However, at short separations, the probabil- 
ity is largely enhanced, a clear evidence of a short-range 
attractive interaction between the two particles. As the 
size of the system is increased, the enhancement peaks 
up, although the interaction contribution to the total 
energy decreases and becomes vanishingly small in the 
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thermodynamic limit MM- 

The spinon-holon interaction has important conse- 
quences on the functional form of the hole spectral den- 
sity, Ah(ui,q). Ah(u},q) does not exhibit Landau's quasi- 
particle resonance 0- It rather shows a sharp square- 
root threshold, followed by a broad branch-cut. Broad 
spectra have also been experimentally seen in ARPES ex- 
periments performed on quasi 1-d samples [[fl| . We prove 
that the divergent square-root threshold is the main ef- 
fect of the attraction between a spinon and a holon. The 
attraction enhances the matrix element for decay of a 
hole into a spinon-holon pair, thus making the hole exci- 
tation fully unstable to decay into a spinon and a holon. 
In the thermodynamic limit, the enhancement turns into 
the sharp square-root threshold, followed by the branch 
cut. 

Let us start our analysis with some basic results from 
the KY-model. At filling 1/2, the ground state of the 
KY-model is the same as the ground state of the HS- 
model. For even N, the corresponding wavefunction is 
given by 
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*gs(zi, . . . , z M ) = y\_( z i - z j) 2 n • 



(2) 
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where M — N/2 and the {j}'s denote the locations of j- 
spins, all the others being j. The corresponding energy is 
given by £; G s = - J(7r 2 /24)(iV+5/7V) A j spinon 

excitation localized at s, at fixed filling, can be thought 
of as a singlet sea where the spin at s is constrained to 
be j. It is described by the wavefunction 
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f s (>l, . . . , Z M ) = Y[(Zj - s)Zj Y[(Zi - Zjf 



(3) 



where now N is odd and M = (N — l)/2. Onc-spinon 
eigenstates of Hky are given by spinon plane waves, 



, zm) 



(4) 



where to = 0, . . . , M and the total crystal momen- 
tum is q% = (ir/2)N ~ {2ir/N)(m + 1/4) (mod 2tt) 
The energy is given by E% = -J{ir 2 /24)(N - l/N) + 
J /2(2tt /N) 2 m(M — to). In terms of q%% the energy with 
repect to the ground state is E(q%) = (J/2)[(tt/2) 2 - 
(C) 2 ] (mod tt). 

One-holon states carry spin and charge 1 with re- 
spect to the ground state. They can be constructed by 
removing an electron from the center of a spinon J7],§],[l4| . 
Unlike in the spinon case, the Brillouin zone for one-holon 
states is not halved juj . The negative-energy part of the 
Brillouin zone is spanned by the propagating one-holon 
wavefunctions 



A I 



A I 



i°{ Z1 , . . . , ZM \h) = h n n( Zj - h)zj - » (5) 



where h is the coordinate of the empty site and < n < 
(N + l)/2. The total crystal momentum of the state in 
Eq.(§ is q%° = (tt/2)N+ (2n/N)(n- 1/4) (mod 2tt). Its 
total energy is given by E%° = -J(n 2 /2A)(N - l/N) + 
J/2(2n/N) 2 n(-(N + l)/2 + n). The energy with respect 
to the ground state is E ho {q^°) = -(J '/2)[(tt /2) 2 ~(q*°) 2 } 
(mod 7r) [p]J^Jl4]|. Here we will not consider positive- 
energy one holon eigenstates |l4|], since they are irrel- 
evant to our analysis. 

Spinons and holons do not lose their identity in a 
many-spinon many-holon state. Hence, to construct one- 
holon one-spinon eigenstates of i?KY we can start with 
states in a mixed representation, where only the spinon 
is localized at s. Let N be even and M = N/2 — 1. The 
state 'J™ is defined as 
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. . . ZM \h) = h n n( Zj ~ S )(zj - h) Zj H(zi ~ Zj ) 2 , 



><j 



(6) 



where 1 < n < M + 2, the spin at s is constrained to be 
I and h is the coordinate of the empty site. States with 
a well-defined crystal momentum can be constructed by 
propagating the spinon and are given by 

* mn (zi,-.-z M \h)=Y,^-^(z 1 ,...z M \h) . (7) 



The total crystal momentum is q = (ir/2)(N -2) + q^ + 
q^° (mod 2n) ; q^ and q^° are the momenta of the sin- 
gle spinon and holon, respectively. One-spinon one-holon 
eigenstates of the KY-Hamiltonian are linear combina- 
tions of the states ^ mn with fixed total momentum 
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The coefficients at are defined by recursion as 
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a k a = 1 



(10) 



and the corresponding eigenvalue is 



i<3 



E mn = E GS + E(q%) + E{ql°) - ^ IC 2 9 " Q| ■ (11) 

The energy of the one-spinon one-holon state relative to 
the ground state is the sum of the energies of an isolated 
spinon and an isolated holon plus an interaction term 
that is vanishingly small in the thermodynamic limit. 
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In order to compute the norm of the states <E> m „, 
we employ the recursive technique introduced in 
§,§. ^ has the generic form ^"(zi, . . . , z\i\h) = 
<3>"(zi, . . . , zm|^)^gs, where $™ is a symmetric polyno- 
mial in z±, . . . , zm, and 'J'gs is the wavefunction intro- 
duced in Eq.(|2|). Define the operator e\(z\, . . . , Zm) — 
Z\ + . . . + Zm- It is straightforward to prove that 



[ J ffKY,e 1 ]K-* GS ]=*GS^(^) 2 
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ds 2 



h- 



(c \ 71 



-M(s + /i) 



(12) 



Using Eq. (|12]) , we find the following relations between 
matrix elements of [Hky, ei] between energy eigenstates, 
calculated with respect to the inner product (f\g) = 
T, Zl ,..., ZM ,hf*( z i>---> z M\h)g(zi,...,z M \h) and valid in 
the case m — n + 1 < 
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($ m -l,n|ei|$mn) _ 

($ m _ 1)n |$ m _ 1)n ) ~ 2(M 



and 



(4> 



m— 1 ,n | ^-1 



/?? 



($m,n\®m, n ) 2(m - |) 

From eqs.(13 |3) one finds by recursion that 

($mn\®mn) _ T[m + \]Y[M - m + 1] 



(*gsI*gs) 7vr[m + i]r[M 



(13) 
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where (#gs|*gs) = N M+1 (2M + 2)!/2 M+1 JlSj . Asm- 
?i + 1 > 0, by following the same steps we find 



T[m+1]T[M 



(*gs|*gs> NT[m + %]T[M- 



(16) 



The state for a localized spinon at s and a localized 
holon at ho, ^ s h is defined as the Fourier-transform of 
'J™ back to coordinate space 



**h = 



M+2 

E v*? 



(17) 



As in [|8|,|9|, we define the real-space coordinate represen- 
tation for a spinon-holon pair, s m h ~ n p mn (s/ho), as 



M+2 M 

* S fc = E E sm V>m"( S / /l o) $ r, 
n— 1 m— 



(18) 



Notice that, from Eqs.(p|,p|Jl^), we find that the relative 
wavefunction, p mn , is a polynomial in the variable s/ho 
if m — n + 1 < 0, and a polynomial in the variable /io/s 
if m — n + 1 > 0. Since $ m „ is an eigenstate of the 
KY-Hamiltonian, we obtain 



(E mn - E GS )($ mn \V sho ) = ($ m „|(# K Y - E GS )\$ sho ) 
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where v = M i£ m — n + 1 < 0, 1/ = otherwise. 

For m — n + 1 < 0, Eq.(|l^) turns into the following 
equation of motion for p mn 
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1 (1- 
whilc, if to — rt 



+ 1 > 0, it becomes 



-p mn (l) = (20) 
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Pmn (^) 



(1)=0. (21) 



Eqs.(^0|,|l]) are Dirac-like first order differential equa- 
tions. They contain a short-range interaction potential 
that diverges at short distances as the first power of the 
separation between spinon and holon. Its effects can be 
analyzed by studying the corresponding exact solutions. 
The solution to Eq. (|20|) is a polynomial in z 
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(22) 



while the solution to Eq. (pl|) is a polynomial in 1/z 
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(23) 



In Fig.|] we plot |p m „(e l9 )| 2 vs. 9. The sharp maxi- 
mum at small spinon-holon separation is the main ef- 
fect of the strong attractive potential in Eqs.( p0| , pT| ). It 
is worth stressing that the interaction between a spinon 
and a holon has, for large N , exactly the same shape as 
the interaction between two spinons H]|J, although the 
differential equation for the two-spinon wavefunction is 
second order, while Eqs.( p2| , p3| ) are first order. 

To rigourously prove that the spinon-holon attraction 
generates the square root singularity followed by a branch 
cut in the hole spectral function, A^(q, uj) [fjj, we now cal- 
culate the contribution to A^q^oo) from one-spinon one- 
holon states. Since i?KY acts on Gutzwiller-projected 
states, matrix elements of Z?ky between states with at 
least a doubly-occupied site are zero. Accordingly, at half 
filling Ah(u>,q) takes contributions only from hole states 
propagating forward. It is given by 
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A h (q,uj) = ^2 
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TrN(X\X)(* GS \* GS ) 
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oj + if] - (E x - Eqs) 



(24) 



where \X) is an exact eigenstate of Hky, Ex is its energy 
and q = 2nk/N. In order to calculate the contribution 
to ^4h(<Z, w) from one-spinon one-holon states, °(q, u>) 
we sum over \X) = |$ mn ). From Eq.(p7|) we obtain 

r M+2 1-2 x 2 /i\ 

ispho/ \ rv 1 I °k-m+lP m i\l) . 
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4-m+i(Pm;) 2 (l) 
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(25) 



This proves that the contribution to the hole spectral 
density from one-spinon one-holon states is fully deter- 
mined by the probability enhancement due to the short- 
range spinon-holon attraction. 

The thermodynamic limit of Eq. ( p5|) is defined as M — > 
co, with m/M and n/M constant. By using the identity 
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J r[i]r[L + i] 
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spinon-holon wavefunction becomes the square-root sin- 
gularity in ho (ui , q) at the threshold energy for cre- 
ation of a spinon-holon pair, followed by the broad branch 
cut. Spinon-holon attraction enhances the matrix ele- 
ment for decay of a hole into a spinon-holon pair. The 
physical consequence is the instability of the hole, which 
is no longer a legitimate excitation of the system. 

In conclusion, we have studied the spinon-holon inter- 
action in an exact solution of the t — J model with 1/r 2 
interaction. A spinon and a holon interact by means of 
a short-range attraction. Although unable to bind the 
two particles, it generates a probability enhancement as 
they are close to each other. It corresponds to an en- 
hancement in the matrix element for decay of a hole into 
a spinon-holon pair, which makes the hole excitation un- 
stable. Hence, our result shows that spinon-holon attrac- 
tion is what makes Landau's Fermi liquid theory break 
down in 1-d strongly correlated electron systems and, 
consequently, makes the quasiparticle resonance disap- 
pear |lr| . 
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and by approximating the gamma functions with Stir- 
ling's formula, the thermodynamic limit of Eq.([25|) can 
be written as an integral over the one-spinon one-holon 
Brillouin zone 



A s v ho (uj,q) = 23m 



v dq } 
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dq s - 



* dq s 



<% ~ gsp + gho) 

7r - q sp J lo - fi + ir\ - E(q sp , q ho ) 



(27) 



(In Eq.(|27|) we have added the chemical potential /i = 
Jtt 2 /4, in order to fix the filling at (N - 1)/2N. qsp 
and q^o are the spinon and holon momenta.) Eq.(|27j) 
has been first worked out in |Q. Af p ho (ui, q) is ^ for 
< q < 7r. In this region of values of q, integration of 
Eq.(^) is straightforward. It provides: 
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From Eq.(p8|) we see that, in the thermodynamic 
limit, the short-distance probability enhancement in the 
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